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ABSTRACT 

We present an analysis of the 0.95—14.5 fim spectral energy distributions of nine field ultracool 
dwarfs with spectral types ranging from LI to T4.5. Effective temperatures, gravities, and condensate 
cloud sedimentation efficiencies are derived by comparing the data to synthetic spectra computed 
from atmospheric models that self-consistently include the formation of condensate clouds. Overall 
the model spectra fit the data well, although the agreement at some wavelengths remains poor due to 
remaining inadequacies in the models. Derived effective temperatures decrease steadily through the 
LI to T4.5 spectral types and we confirm that the effective temperatures of ultracool dwarfs at the 
L/T transition are nearly constant, decreasing by only ~200 K from spectral types L7.5 to T4.5. The 
condensate cloud properties vary significantly among the L dwarfs in our sample, ranging from very 
thick clouds to relatively thin clouds with no particular trend with spectral type. The two objects 
in our sample with very red J — Ks colors are, however, best fitted with synthetic spectra that have 
thick clouds which hints at a possible correlation between the near-infrared colors of L dwarfs and 
the condensate cloud properties. The fits to the two T dwarfs in our sample (T2 and T4.5) also 
suggest that the clouds become thinner in this spectral class, in agreement with previous studies. 
Restricting the fits to narrower wavelength ranges (i.e., individual photometric bands) almost always 
yields excellent agreement between the data and models. In some cases, the models even reproduce 
the detailed structure of the weak II2O and CH4 absorption features. Limitations in our knowledge 
of the opacities of key absorbers such as FeH, VO, and CH4 at certain wavelengths remain obvious, 
however. The effective temperatures obtained by fitting the narrower wavelength ranges can show a 
large scatter compared to the values derived by fitting the full spectral energy distributions; deviations 
are typically ~200 K and in the worst cases, up to 700 K. In spite of the uncertainties remaining in 
the models, the approach presented here works rather well and should be extended to a much larger 
sample. 

Subject headings: infrared: stars — stars: low-mass, brown dwarfs stars: fundamental — radiative 
transfer 



1. INTRODUCTION 

The emergent spectra of very low-mass (VLM) stars 
and brown dwarfs, collectively known as ultracool dwarfs, 
are shaped in large part by gas and condensation chem- 
istry The low temperature (500 K < T < 3000 K), 
high pressure (1 bar ^ P ^ 10 bar) conditions typi- 
cal of ultracool dwarf atmospheres favor the formation 
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of molecules (such as CO, CH4, N2, NH3, and H2O) 
and, at T < 2400 K condensates from the refractory 
elements (Ti, V, Ca, Al, Fe, Si, and Mg). Condensa- 
tion removes gaseous opacity sources (such as TiO, VO, 
FeH, CrH) from the atmosphere and alters the atmo- 
spheric chemistry as the condensates gravitationally set- 
tle and form clouds that deplete the upper layers of the 
atmosphe re of both the condensate s and their constituent 
eleme nts (iFedev fc LoddersI 11994 119961 : iLoddersI I1999L 
l2002t iBurrows fc Shard Il999l )^ Condensates also con- 
tribute their own opacity, which can be substantial for 
the more abundant Fe-, Mg-, and Si-bearing species such 
as forsterite ( Mg2Si04), enstatite ( MgSiO.s), a nd sohd or 
liquid Fe (e.g. jMarlevet al.ll2002l:lTsuiill2(302[ ). 

This rich chemistry gives rise to complex emergent 
spectra that are distinctly non-Planckian. The molecular 
line opacities dominate over all of the continuum opac- 
ity sources (e.g., H~) t hat are important in the a tmo- 
spheres of hotter stars (|Allard fc Hauschildtl Il995f ) and 
as a result the emergent spectra of ultracool dwarfs ex- 
hibit strong molecular absorption bands of TiO, H2O, 
and CH4. Indeed, it is the strength of these and other 
molecular band s that are used to define the L and T 
spectr al classes ([Kirkpatrick et al.]|1999l : iBurgasser et al.l 
[2006al) . L dwarf (2400 K < Tes < 1400 K) spectra ex- 
hibit weakening TiO and VO bands at red optical wave- 
lengths (0.6—1.0 /im) due to the formation of Ti-bearing 
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condensates such as perovskite (CaTiOs). The loss of the 
gaseous TiO and VO opacity results in a more transpar- 
ent atmosphere and leads to the emergence of monatomic 
alkali lines (Na, K, Rb, Cs). Most prominent are the res- 
onance Na I (5890/5896 A) and K I (7665/7699 A) dou- 
blets, which become increasingly broad with later spec- 
tral type. At near-infrared wavelengths (1—2.5 ^m) the L 
dwarfs also become progressively redder with later spec- 
tral types due to the additional condensate opacity. As 
CH4 becomes the dominant carbon-bearing gas in the up- 
per, coolest layers of the atmosphere (CO/CH4 < 1 for 
T < 1100 K at P=l bar), CH4 absorption bands emerge 
at near-infrared wavelengths which signals the transition 
between the L and T (1400 K < Teg < 700 K) spectral 
classes. The K I and Na I resonance l ines are present 
in the red optical spectra of T dwarfs (jBurgasser et al.1 
l2003a^ because the high temperature condensate anor- 
thite (CaAl2Si208), a feldspar, has gravitationally set- 
tled out of the cooler atmospheric regions. In the ab- 
sence of such settling, Na and K would condense into 
a feldspar in the form of albite (NaAlSiaOs) and ortho- 
clase (KAlSi.^Os) ([Burrows fc Shardll999l : lLodderslfl999l : 
iMarlev et aLil200a K 

The observed spectra and photometry of ultracool 
dwarfs can therefore serve as touchstones against which 
we can test our understanding of atmospheric physics, 
chemistry, dynamics, and cloud processes as embodied 
in our atmospheric models. Shortly after the discov- 
ery of the first bona fide brown dwarf, Gliese 229B 
(jNakaiim a et al.l Il995f ). a number of publications ap- 
eared comparing the observed spectrum to model s 



■Allard et al.lll996HMarlev et al.lll996l : lTsuii et al.lll996D . 
iBurrows et al.l ()2001l ) presented a review of the substan- 
tial subsequent work devoted to understanding the prop- 
erties of L and T dwarfs. Although the comparison of 
models to the spectra of individual objec ts such as Gl 
229B (Saumon ct al. 2003) and Gl 570D (iGeballe et all 
I2OOII : ISaumon et al.ii2006l ) has played an important role 
in the progression of the field, relatively little work has 
been done on comparing entire classes of models to uni- 
form datasets. Such comparisons are important as they 
can reveal the global properties of a set of objects and 
highlight the systematic strengths and weaknesses of the 
models. 

In this work, we present such a comparison to one 
particular set of atmospheric models (Marley et al. 2002, 
2007, in preparation). We aim to understand what the 
modeling has accomplished to date and to help direct 
future models and observations along the pathways that 
seem most promising. The paper is organized as fol- 
lows. Sections §2 and §3 present the data and model 
atmospheres used in the analysis, respectively. In §4, we 
present the results of the comparisons between the data 
and models and in fj5l we discuss two systematic uncer- 
tainties that affect our analysis. 

2. THE DATA 

There are currently eleven L and T dwarfs that 
have published spectra with nearly complete wave l ength 
coverage from 0.65 to 14.5 /xm (jGushing et al.1 l2006t 
ISaumon et all 120071). We sele cted the 5 L and 2 T 
dwarfs from the iGushing et al.l sample with the high- 
est S/N spectra that are not known to be multiple sys- 
tems. We also included two additional L dwarfs, 2MASS 



J1506514-I-1321060 and DENIS J02550 3.3-470049.0, 
that were not part of the IGushing et al.l sample. The 
properties of the dwarfs in our sample, including opti- 
cal and infrared s pectral types and Tw o Micron All Sky 
Survey (2MASS; 'skru tskie et al.ll2006l ) photometry, are 
given in Table [TJ Although both optical and infrared 
spectral types are listed in Table [H we hereafter us e op- 
tical types for the L dwarfs (jKirkp atrick et al."1999') and 
infrared types for the T dwarfs (Burgasser et al. 2006a). 
In addition, we abbreviate the numerical portions of the 
2MA SS, Sloan Digital Sky Survey (SDSS; York ct all 
|2000( ). and Deep Near Infrared S urvey of the Southern 
Sky (DENIS : lEpchtein et anil997| ) target designations as 
hhmmibddmm, where the sufhx is the sexigesimal Right 
Ascension (hours and minutes) and declination (degrees 
and arcminutes) at J2000 equinox. 

The 0.65—14.5 /im spectra used in this work were con- 
structed from the following sources: red-optical spectra 
were obtained wit h the Low Resol ution Imaging Spec- 
trograph (LRIS; lOke et al.l Il995| ) at the Keck tele- 
scope; near-infrared spectra were obtained with SpeX 
(jRavner et al.ll2003[) at the NASA IRTF and the Infrared 
Camera and Spectrograph (IRCS; iKobavashi et al.l 
I2OOOD on Subaru; and mid-infrared (5.5—14.5 fim) spec- 
tra were ob tained with the Infrared Spectrograph (IRS; 
iHouck et al.li2004i) o n board the Spitzer Space Telescope 
(jWerner et al.ll20M ). The only wavelengths not covered 
by these spectra are those containing the strong telluric 
absorption bands at 1.9 and 2.6 /Ltm and the 4.5—5.5 fim 
range. The latter wavelength range is exceedingly diffi- 
cult to observe from the ground and is not covered by 
the IRS. Table |2| lists both the references and resolv- 
ing powers of the various spectra used in this work. For 
each dwarf, the different spectra were absolutely flux cal- 
ibrated using ground-based photometry and then com- 
bined into a single 0.65—14.5 /ini spectrum as described 
in IGushing et al.1 (|2005, 2006). 

Figure E shows the 2MASS J-H,H-K„ and J-Ks 
colors as a function of spectral type for all of the L and T 
dwarfs listed in the Dwarf Archives^^ web site as of 2007 
June 07. The dwarfs in our sample are shown as red 
points. Overall, our sample is representative of dwarfs 
with similar spectral types. However, 2MASS 2224-0158 
(L4.5) and 2MASS 0825+2115 (L7.5) have relatively red 
J — Kg colors given their spectral types and 2MASS 
J1507-1627 (L5) has a relatively blue J - K^ color. We 
discuss these dwarfs in more detail in §4.61 

3. THE MODELS 

The L and T dwarf models presented here have 
been produced in the same way as those in our ear- 
lier work (e.g.. iMarlev et al.l [19961: IBurrows et al.lll997l : 
IMarlev etal1l2002l : ISaumon et al.ll2006[ ). We are prepar- 
ing a complete description of the models for a fu- 
ture publication. Briefly, we computed a grid of 
one-dimensional, plane-parallel, hydrostatic, non-gray, 
radiative-convective equilibrium atmosphere models us- 
ing our standard modeling code. A one-dimensional, 
plane parallel approximation is appropriate since the at- 
mospheric scale height is always much smaller than the 
ultracool dwarf radius. We do not use mixing length 
theory, but rather set the local lapse rate in the con- 

^^ |http : //Dwarf Archives ■ org | 
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vective regions to the adiabatic lapse rate. This is an 
excellent approximation in the dense, warm atmospheres 
of brown dwarfs and giant planets where convection is a 
very efficient energy transport mechanism. The supera- 
diabaticity is very small and the atmospheric pr ofile is 
insensitive to the choice of mixing length (Baraffe et al.l 
fl997t lAHaTd et aLlHool . 

For the chemical equilibrium calculations we use the 
elemental abundance data of Eoddcrs (2003) an d corn - 
pute composition s following Fcgley & Loddcrs (1994'), 
iLodders fc FeglevI (2002 ) . and Lodders (1999. 2002). The 
chemical equilibrium abundances are computed for lo- 
cal thermodynamic equilibrium conditions. In stellar at- 
mospheres, departures from thermochemical equilibrium 
can arise from interactions of atoms and molecules with 
the non-thermal radiation field while in planetary atmo- 
spheres, non-equilibrium conditions can arise from ver- 
tical transport (by convection or eddy diffusion) on a 
time scale shorter than that of the chemical reactions in- 
volved. In ultracool dwarf atmospheres, the first effect is 
negligible compared to the second. In the present work, 
we neglect both effects. The opacity database includes 
the molecular lines of H2O, CH4, CO, NH3, H2S, PH3, 
TiO, VO, CrH, FeH, CO2, HCN, C2H2, C2H4, and C2H6, 
the atomic lines of the alkali metals (Li, Na, K, Rb and 
Cs), and continuum opacity from H2 collision induced 
absorption (CIA), Rayleigh scattering from H2, H and 
He, bound- free opacity from H, H~, and H^, and free- 
free opacity from He, H^, and H^. The molecular line 
lists are continually updated and while they have well- 
known shortcomings, they represent the state-of-the-art. 
An explanation of the treatment of line broadening, and 
further details regarding the chemical equilibriu m cal- 
culation are presented in iFreedman et alT (|2007f ). Note 
that we do not include a turbulent velocity in the cal- 
culation of line shapes as test calculations have shown 
this has a negligible effect for the temperatures and pres- 
sures consi dered here. An error in the TiO line opacity 
database (Schwenkelll998[ ) employe d in our recent gener - 
ations of models (including th ose in | Marlev et al. (2002f) , 
iBurgasser eFall (|2002b.) . and lKnappet all (|2004.) ) that 
resulted in too-strong TiO absorption has been corrected. 

For the model atmosphere calculation we use the k- 
coeffi cient method to describe the opacities (|Goodv et al.l 
|1989() . Within each of 180 spectral bins we sum the 
opacity arising from each of the several hundred mil- 
lion atomic and molecular lines in our database. Our 
radiative transfer follows the source function technique 
(|Toon et al.l Il989f ) , allowing inclusion of arbitrary Mie 
scattering particles in the layer opacities. This ap- 
proach is exact in the limit of single scattering and treats 
multiple scattering by assuming that the scattered ra- 
diation field can be approximated with two streams. 
After computing a pressure-temperature (P, T) profile 
with the fc-coefficient method for a given set of param- 
eters, we compute a high resolution spectrum by solv- 
ing the monochromatic radiative transfer equation for a 
larger number of frequency points (typically ~ 190000 
in the range of 0.4 to 50 /zm), using the exact same 
molecular line lists, continuum opacity sources, chemi- 
cal composition, and line profiles as those used for the 
computation of the fc-coefficients. The resulting spec- 
tra can be smoothed or binned for comparison with 



data. Previous applications of our atmospheres code 
include the generation oi P/T profil es and spectra for 
Titan (jMcKav et a ll '1989!) , Uranus (iMarlev fc McKavl 
I1999D . hot Jupiters (Fortn ev et ap l2005l I2006D and 
brown dwarfs (jMarlcv et al. 1996; B urrows et al.l fl997t 
IMarlev et all[2002 : Saumon et al. 20061 ). 

As we converge to a model temperature-pressure pro- 
file, we compute the conden sate size and vertical profiles 
following the prescription of lAckerman fc MarlevI ()2001l ) 
for various values of the cloud sedimentation efficiency 
parameter, fsed- Larger values of /sod imply larger par- 
ticle sizes and thus greater sedimentation efficiency. The 
width of the lognormal particle size distribution is set at 
cr = 2 and the eddy diffusion coefficient above the con- 
vective region is set to K^z — 10^ cm^ sec^^. This latter 
value is consistent with that found from studies of the 
effect of no n-equilibrium chemistry in L and T dwa rf at- 
mospheres ([SaumoZiEaDllQQi EeggittiranilOO^) . For 
any P/T profile there is a single description of the varia- 
tion in particle sizes and number densities with altitude 
for each of the condensed species. As the atmospheric 
structure is iterated towards the equilibrium solution the 
cloud description is also continuously updated. The final 
result is a single, self-consistent thermal and cloud pro- 
file with the specified gravity, effective temperature, and 
cloud sedimentation efficiency. We term these models 
'cloudy'. We also consider models in which the conden- 
sate formation is included in the chemical equilibrium 
calculation, but the condensate opacity is neglected. We 
term these models 'clear' and denote them as /sod=nc, 
since the models include 'no cloud' opacity. 

In this first study, for the sake of simplicity, we consider 
a model parameter space consisting of Teff , g, and /sod at 
solar metallicity A more extensive study including vari- 
ations in metallicity and the other cloud parameters, as 
well as departures from chemical equilibr ium caused by 
vertical mixing in the atmosphere (e.g. ISaumon et al.l 
(|2006D ). will be the subject of future papers. Evolu- 
tionary models predict that field L and T dwarfs with 
ages greater than a few hundred million years have 700 
K < Toff < 2400 K and 5.0 < logg (cm s'^) < 5.5 
(jBurrows et al.l 119971 ). We consider the entire range of 
Toff in steps of 100 K, and values of logg of 4.5, 5.0, 
and 5.5. In previous work found that /sod ~ 3 best re- 
produces global properties of Jupiter's ammonia cloud 
(Ackerman & Marlev 2001) and the J — K colors of L 
dwarfs (e.g., iKnapp et al.l 120041 ). However the excess 
TiO opacity in our previous models, particularly in the 
J band, may have led us to overestimate the fsed value 
for L dwarfs. Accordingly here we consider models with 
/sed=l, 2, 3, 4, and nc. The final grid contains 270 mod- 
els. 

Dynamical processes in the atmospheres of individ- 
ual objects ultimately control the growth and sedi- 
mentation of condensates and the optical properties of 
the v arious cloud decks e xpected in these atmospheres 
(e.g., LH clhng et al.' 2001; Woitkc & Helling 2003, 2001; 
iHelling et al. 2004 : Helling fc Woitkc 2006 ) . Although 
our modeling parameter, fsed aims to self-consistently 
capture these processes in a tractable way, the global 
properties of each individual object must ultimately de- 
fine the true properties of the clouds. However we do not 
yet understand how differences in metallicity, rotation 
rate, effective temperature, and gravity might manifest 
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themselves as different values of /sod- Therefore while 
we treat this parameter on equal footing as g and Teg 
for fitting purposes, we recognize that ultimately, more 
complete models would find that /sed is determined by 
the global physical properties of individual objects. 

For each combination of gravity, effective tempera- 
ture, and sedimentation efficiency, we compute a self- 
consistent atmospheric temperature-pressure profile. Us- 
ing this profile and employing the same atmospheric 
chemistry and opacities, we compute a high resolution 
synthetic spectrum. This spectrum is then smoothed to 
the resolving power of the data using a Gaussian kernel 
and resampled onto the same wavelength grid. Figures [2] 
and [3] show a sequence of synthetic spectra, smoothed to 
i?=500 for clarity, with atmospheric parameters typical 
of an L dwarf (Toff =1800 K, log5=5.0, /scd=2) and a 
T dwarf (reff=1100 K, log5=5.0, /sod=nc). The three 
panels in each figure illustrate the variations in spec- 
tral morphology when just one of the three parameters 
is changed. 

The upper panels show the changes in spectral mor- 
phology due to variations in T^s at a fixed g and /sed- 
With decreasing Tefi, the L dwarfs (Fig. [2]) become 
redder in the near-infrared and the depths of the H2O 
bands become deeper. Note also that the Q branch of 
the v^ fundamental band of CH4 at 3.3 /im is present 
at Toff =1700 K. The depths of the H2O, CH4, and NH3 
bands in the spectra of the T dwarfs (Fig. [3]) become 
deeper with decreasing Tcff but in contrast to the L 
dwarfs, the T dwarfs become slightly bluer in the near- 
infrared. 

The middle panels show the changes in spectral mor- 
phology due to variations in g at a fixed Tgff and /sed- 
In the L dwarfs the primary effect of increasing g is 
to weaken the near- infrared H2O absorption bands. At 
these temperatures, the molecular abundance of H2O and 
GH4 depend weakly on the atmospheric pressure and are 
thus not sensitive to changes in gravity. The change in 
the H2O band strengths is instead a result of the fact 
that, at a given pressure level in the atmosphere, the 
higher gravity models are cooler and as a result have 
thicker condensate clouds. The clouds produce a grayer 
emergent spectrum and are therefore responsible for the 
decreased depth of the H2O bands. In contrast to the 
L dwarfs, the H2O and GH4 band depths in the spectra 
of T dwarfs are only marginally affected by changes in 
g, except in the K band. The peak flux level in the K 
band diminishes with increasing g due to an increase in 
the G IA of H2 - This broad , featureless absorption band 
(e.g., iBorysow et al.l [19971 ) centered at 2.4 ^m is very 
sensitive to atmospheric pressure, and thus g, since it in- 
volves the collisions of particles (k^, oc "n? g, where Ugas 
is the number density of the gas). 

Finally, the lower panels show the changes in spectral 
morphology due to variations in the cloud sedimentation 
parameter /sod at a fixed Tcff and g. Ghanges in /sod 
produce the most dramatic variations in spectral mor- 
phology (particularly in the 1—6 /im region) underscor- 
ing the importance of correctly modeling the formation 
and subsequent sedimentation of condensates. The ad- 
ditional condensate opacity results in emergent spectra 
that are redder at near-infrared wavelengths than those 
from the cloudless models (/sod=nc). The impact of the 



clouds on the emergent spectra of ultracool dwarfs be- 
comes more dramatic with lower /sed values because the 
clouds are thicker and thus contribute more opacity to 
the atmosphere. 

4. ATMOSPHERIC PARAMETERS OF L AND T DWARFS 
4.1. Fitting Technique 

The ID model atmospheres predict the emergent flux 
density Ti, ( Jv = 27r J^ /^/x d/i) as a function of wave- 
length (i.e., a synthetic spectrum) at the surface of a star. 
The synthetic spectra must be multiplied by {R/d)^, 
where R is the stellar radius and d is the stellar dis- 
tance, in order to obtain fluxes at Earth and directly 
compare them to observations. Although the parallaxes 
of nearly 80 L and T dwarfs have been measured (e.g , 
Dahn e t al.' '2002; ' Tinnev et all [200I IVrba et all [200I 
Gosta et al. 2006) , the radii of field L and T dwarfs are 
currently unknown. Therefore we can compare only the 
relative shapes of the modeled spectral energy distribu- 
tions (SEDs) to the data. 

For each model k, which hereafter is denoted by the 
triplet [Toff/log 5//sod], we compute a goodness-of-fit 
statistic Gfc defined by^^. 



G^-E- 



i=l 



fi — Ck^kA 



(1) 



where n is the number of data pixels, Wi is the weight 
for the ith wavelength, fi and !Fk,i are the flux densities 
of the data and model k, respectively, ai are the errors 
in the observed flux densities, and Ck is the unknown 
multiplicative constant equal to {R/dY . For each model 
fc, the constant Ck is determined by minimizing Gk with 
respect to Ck and is given by. 



Ck 






(2) 



We select the best fitting synthetic spectrum by locating 
the global minimum of the G values for all of the model 
spectra in our grid. 

We have no a 'priori reason to favor one wavelength 
range over another and therefore, in principle, all wave- 
lengths points could receive equal weight in our fits. 
However given the large variation in the wavelength sam- 
pling of our data (e.g, the near-infrared spectra cover 
only ^2 /im but have ~20 times more wavelength points 
than the Spitzer IRS spectra which cover ~10 /xm), 
G would be heavily biased towards near-infrared wave- 
lengths if Wi was set to unity for all i. We therefore 
chose to weight each pixel by its width in microns {wi — 
AXi). In i i5.2l we discuss the systematic effects that our 
choice of weights has on our results. Since the line lists of 
GH4 are incomplete in the 7J-band (jBurrows et al.ll2001l : 
iFreedmanet al.ll2007f) . and the E '^U - A '^U system 
of FeH at 1.6 /im which is present in the spectra of the 



We also tested the usefulness of the statistic Gt 



En 
i=l-^i 



fi 



since most published spectra do not have 



er rors associated with the flux density values. However as noted 
bv lTakedal 119951) , such a statistic gives more weight to wavelengths 
with low flux values. Typically ultracool dwarf spectra have very 
low S/N at such wavelengths and thus care must be taken when 
using such a statistic. 
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L dwarfs (|Wallace fc Hinkld[200ll : iCushing et aT][200l 
lacks a line list, we also set the weights of the pixels from 
1.58 to 1.75 /im (roughly the center of the H band) to 
zero. 

As an example of our procedure, Figure [4] shows the G 
values of all the synthetic spectra in our grid as a function 
of Toff computed for the spectrum of 2MASS 1507-1627 
(L5). In this case, the best fitting synthetic spectrum 
has the parameters [1700/4.5/2]. However there are ad- 
ditional synthetic spectra that yield similar minimum G 
values. How then can we determine the quality of a given 
fit? Although G is mathematically similar to the x^ 
statistic, it does not follow a x^ distribution^'^ and as 
a result, we cannot employ standard hypothesis testing 
techniques to assess the quality of the fits. 

We therefore performed a Monte Carlo simulation to 
determine the range of synthetic spectra that fit the data 
given the observational errors. For each observed spec- 
trum, we generated 2000 simulated data sets wherein the 
flux at each wavelength point was randomly drawn from 
a Gaussian distribution centered on the observed flux, 
and with a width given by the observed variance in the 
individual pixel. Since each 0.65—14.5 fim spectrum was 
constructed by combining together spectra that were ab- 
solutely flux calibrated independently, the error in any 
given flux density value also includes a correlated com- 
ponent that arises from the fact that during the absolute 
flux calibration process, all of the flux density values over 
a range of wavelengths (e.g., 1—2.5 /xm) are scaled by a 
constant to adjust the overall flux level to match broad- 
band photometry. We therefore scale the simulated flux 
density values by IQ"*^-^^^ where K is randomly drawn 
from a Gaussian distribution centered on zero, and with 
a width given by the variance of the photometry (in mag- 
nitudes). 

A simulated data set s^ is constructed from m spectra, 
each having its own a^f^. For a given spectrum j, the 
simulated dataset Si is given by, 



= AA(/„af)10-°-4x^(o>-?'.,.), 



(3) 



where Af{fi, cr^) represents a value randomly drawn from 
a Gaussian distribution with mean /x and variance cr^, 
and a\ ■ is the variance of the photometry used to ab- 
solutely flux calibrate the spectrum j. We then replace 
fi with Si in Equations 1 and 2 and recompute the best 
fitting synthetic spectrum. The quality of the fit between 
the data and a given synthetic spectrum is given by /mCj 
the fraction of the 2000 simulations in which the given 
synthetic spectrum is identified as the best fitting spec- 
trum. In the case described above, the synthetic spec- 
trum with the parameters [1700/4.5/2] has /mc=1-000 
and thus is the best representation of the data in our 
grid of synthetic spectra. Values of /mc much less than 
unity imply that more than one synthetic spectrum in 
the grid fits the data well. 

The best fitting model parameters are given on the 
model grid points without interpolation. Nominally, this 

^^ The x^ distribution with u degrees of freedom is defined to 
be the sum of u squared standard normal deviates or x^ = y~! z? , 
where Zi = {xi — ^)/ai. Since the models are known to have 
systematic errors, the standardized residuals Zi = (fi — C}^T}^^i)/ai 
are not standard normal deviates and thus G does not follow a x^ 
distribution. 



would give an internal uncertainty of half of the grid 
spacing, or ±50 K in Tcff, ±0.25 in log^, and ±0.5 in 
/sod- In practice, the best fitting models nearly all have 
/mc^ 1; which shows that the internal uncertainties on 
the fitting parameters are even smaller than those nom- 
inal values. As we will see in [J5l our fits of the ultracool 
dwarf SEDs are affected by systematic effects that dom- 
inate the internal uncertainties and that are difficult to 
quantify. In addition, changes in the models (e.g., the in- 
clusion of [Fe/H] as an additional parameter) or the use 
of a different goodness-of-fit statistic could yield different 
atmospheric parameters than we find here. An accurate 
determination of physical parameters of ultracool dwarfs 
will therefore require the development of more complete 
model atmospheres. Nonetheless, trends in the fitted pa- 
rameters as revealed by this work should improve our 
understanding of the physics of ultracool dwarfs. 

4.2. Fits of the Entire SED 

We fit 0.95—14.5 /im spectrum of each dwarf as de- 
scribed in the previous section. Although the data ex- 
tend blueward to 0.65 /J,m, the analysis was limited to 
wavelengths greater than 0.95 /im because the pressure- 
broadened wing of the resonant K I doublet (7665, 7699 
A) that exten ds to ^0.9 /xm is exceedingly difficult to 
model (Burro ws fc Volobuvev|[2003 ). Table [3] lists the 
best fitting model parameters (Teff, log 5, /sed) for each 
dwarf in our sample along with the fraction of the Monte 
Carlo simulations (/mc) hi which the given triplet was 
identified as the best fitting model. All models with /mc 
> 0.1 are listed. A summary of the derived parameters 
is shown in Figure [5l Seven of the nine dwarfs in our 
sample have best fitting synthetic spectra with /mc "^1 
indicating that they are the only reasonable model spec- 
tra in our grid for these dwarfs. 2MASS 1439+1929 (LI) 
has two model spectra that differ by 1 in /sod with nearly 
equal /mc values while SDSS 1254-0122 (T2) has two 
models with /mc > 0.1 that differ by 0.5 dex in log 5 and 
1 in /scd- 

The derived Toff values decrease steadily with later 
spectral type but the difference in Toff between 2MASS 
0825+2115 (L7.5) and 2MASS 0559-1404 (T4.5) is only 
200 K which confirms that the transition between the 
L and T dwarfs occurs over a nar row range in Tpff 
as has been noted previously (e.g., 'Kirkoatrick et alj 
1999; Burgasser c t al. 2002^ iGolimo wski et al. 2004). 
Golimowski et al.l (J2004l ) derive the effective tempera- 
tures of ~40 L and T dwarfs using their observed 
bolometric l umino s ities and e v olutio nary models (see 
iDahn et al.l (|2002D . IVrba et all (|2004f) . and iBasri et all 
(j2000t ) for alternative Teff scales for L and T dwarfs). 
Seven of these dwarfs a re included in our sample so we 
also show in Table [3] the IGolimowski et all Teff values for 
ages of 0.1, 3, and 10 Gyr. The agreement between the 
two sets of values is excellent as all of our derived values 
fall within their 0.1 to 10 Gyr ranges. 

The comparison of our derived Toff values with those of 
IGolimowski et aLl implies that if they are single objects, 
seven of nine dwarfs in our sample are much younger than 
3 Gyr. Our evolutionary calculations (Marley et al. 2007, 
in preparation) also indicate that if we take the derived 
Toff and g values at face value, five of the L dwarfs would 
be younger than 0.2 Gyr. This is primarily a consequence 
of the low gravities obtained from fitting the SEDs. It 
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seems unlikely that our sample would be biased towards 
such young ages given that the mean age of the L and T 
dwarf field popu lation is several Gyr (|Dahn et al.ll2002l : 
lAllen et al.ll2005f ). This suggests that fitting the SEDs of 
ultracool dwarfs is a poor method to determine gravities. 
Fortunately, we can obtain the gravity using an alternate 
method. 

When the parallax of an object is known, its radius 
can be determined using the fitted normalization con- 
stant (Eq. [2) since Cfc=(i?/d)^. With ultracool dwarf 
evolution sequences, {Tcff,R) values can be transformed 
uniquely into (Teg,g) values. The gravities obtained with 
this alternate method are given in Table [3] and are shown 
in blue in the upper panel of FigurelHl The scatter in log g 
is reduced and nearly all gravities are above log g— 5.0. 
Our evolution sequences show that all 9 objects now have 
ages above 0.3 Gyr. This alternate method only weakly 
depends on the relatively subtle variations in spectral 
morphology arising from changes in g (Fig. [5] and [3]) be- 
cause it relies primarily on fitting the absolute flux level 
of the SED. The reduced scatter in log g and the older age 
distribution of our sample indicate that the method pro- 
vides gravity estimates that are more reliable than those 
obtained by direct spectral fitting. Hereafter, we do not 
compute new model atmospheres with these new gravi- 
ties but rather continue to compare the data to models 
in our grid with log 3=4.5, 5.0, and 5.5. 

We also note that the theoretical evolution of ultracool 
dwarfs with cloudy atmospheres is limited to log(7<5.38, 
although the evolution with cloudless atmospheres can 
reach slightly higher gravities of logg^5.48. Because the 
mass-Tcff-radius relation of brown dwarfs is well estab- 
lished on theoretical grounds, it provides a firm boundary 
to physically plausible combinations of Toff and g. We 
find that the logg values of three of our objects (derived 
using the alternate method) exceed the limit for ultracool 
dwarfs with cloudy atmospheres, which is an indication 
that their fitted Teg values are overestimated by ~10— 20 
K. This result can be understood by noting that the grav- 
ity of brown dwarfs is a very sensitive functi on of Teg for 
a fixe d Lboi (see, for example, Fig. 1 of .Saumon et al.l 
()2006l )). This is a consequence of the mass-radius rela- 
tion for degenerate stars where R ~ M~^^^. Since Lboi 
~ R-^T^ff, it follows that at constant Lboi, g ~ T;^^- From 
this approximate relation, a 7% variation in Tag changes 
the gravity by 0.3 dex. 

Two objects in our sample call for a more detailed dis- 
cussion. SDSS 1254-0122 (T2) and especially 2MASS 
0559—1404 (T4.5) sta nd out as being over luminous for 
their spectral types (jBurgasseil l2007| ) . It has there- 
fore been sugg ested that the two objects may be unre- 
solved binaries (iGolimow ski et a l.ll2004l ; IVrba et al.ll2004l : 
iLiu et al.ll2006l : lBurgasser..2007.) . although neither dwarf 



has been resolved into a bin a ry in h igh-resolution imag- 
ing (Burgasscr ct al.' '2003b|, l2006b[ ). We find, as did 
[Golimowski et al. (2004), that the two objects share es- 
sentially the same T^fF, alt hough we derive Toff=1200 
K while iGolimowski et all find Teff=1475 K (assum- 
ing ages of 3 Gyr). If SDSS 1254-0122 and 2MASS 
0559—1404 are equal magnitude binaries, then the iboi 
of each componen t is reduced by a factor of two and the 
IGolimowski et aLl T^s value becomes ^1200 K which is 
in much better agreement with our value. Given these 



new parameters, our evolutionary sequences provide the 
corresponding properties for the individual components. 
To satisfy both the constraints of component luminosity 
and the maximum allowed gravity, Toff has to be reduced 
slightly, but remain within the uncertainties of our spec- 
tral fits. The properties of the individual components are 
given in Table [H 

Whether SDSS 1254-0122 and 2MASS 0559-1404 are 
singles or equal mass binaries, they have essentially the 
same effective temperature. The difference in the appear- 
ance of their spectra (T2 versus T4.5), can be explained 
as primarily due to a difference in their cloud proper- 
ties, as suggested by our fits that indicate the clouds 
are thinner (larger fsed) in 2MASS 0559-1404 than in 
SDSS 1254-0122. Another possibility would be for SDSS 
1254—0122 to be a single brown dwarf (or a strongly 
unequal-mass binary) and for 2MASS 0559-1404 to be 
an equal mass binary. This would make the latter's indi- 
vidual components less luminous than SDSS 1254—0122 
(or its primary) and is consistent with their position in 
the Mboi vs spectral type diagram dBurgasser 200 7). 

A fair amount of scatter is found in the fsed values 
of the L dwarfs, with an indication that fsed increases 
toward the later T spectral class, as would be expected 
if clouds play a decreasing role in shaping the SED of 
T dwarfs (e.g., Marley et al. 2002; Burrows ct al. 2006). 
Among the L dwarfs in our sample, /sod ranges from 1 
to 3, implying that the cloud opacity varies appreciably 
among L dwarfs (bottom panel of Fig. 2) with little, 
if any, dependence on the spectral type. It would be 
particularly interesting to establish the distribution of 
fsed as a function of spectral type with a larger sample. 

Figure [6] shows the 0.95—14.5 /im spectra of 
2MASS 1439-1-1929 (LI), 2MASS 0036-fl821 (L3.5), 
2MASS 1507-1627 (L5), DENIS 0255-4700 (L8), SDSS 
J1254-0122 (T2), and 2MASS 0559-1404 (T4.5) along 
with the best fitting model spectra. The grey regions 
indicate wavelength ranges that were excluded from the 
fits. The models fit the data reasonably well at the ear- 
liest and latest spectral types, but do a poorer job of 
fitting the data in the mid- to late-type L dwarfs and 
the early-type T dwarf. This range of spectral types cor- 
responds to the L/T transition where condensate clouds 
have their largest impact on the SED of ultracool dwarfs. 
The relatively poor fits at the L/T transition most likely 
reflects the limitations of our simple cloud model. 

In particular, the best fitting model spectra match 
both the peaks and depths of the H2O bands centered 
at 1.4, 1.9, and 2.6 /um well. The models do, however, 
under estimate the ii'-band flux density levels of 2MASS 
1439+1929 (LI) and 2MASS 0036-f 1821 (L3.5). In ad- 
dition, the models also over estimate the depth of the L- 
band CH4 absorption bands of both DENIS 0255-4700 
(L8) and SDSS 1254-0122 (T2) as well as the 5.5-14.5 
/im flux levels of DENIS 0255-4700 (L8). We defer dis- 
cussion of the possible causes of these mismatches to ij4.3l 
Finally, the poor match between the model spectrum and 
data of 2MASS 0559-1404 in the H band is a resuh of 
the woe fully incomplete line list of GII4 at these w ave- 
lengths (jBurrows et al.ll200lHlYeedman et al.ll2007[ ). 

4.3. Description of Individual Bands 

In the following sections, we describe how well the 
best fitting synthetic spectra match the data over nar- 
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rower wavelength ranges. Figures [7] to [12] show the 
same data and model spectra presented in Figure [6] but 
over the wavelengths covered by the photometric band- 
passes Y, J, H, K, and L and the Short-Low mod- 
ule of the IRS on board Spitzer (hereafter IRS/SL). 
Although the most prominent features in the spectra 
of the L and T dwarfs are briefly described in each 
section, more det ailed de scriptions can be found else- 
where (iKlr kpatric k et al.1 \T9 99: Burgas ser et al.1 l2003al : 
iMcLean et al.ll2003l: ICushing et al.ll200al2006f) . 

4.3.1. The Y Band 

Figure [7] shows the spectra and best fitting models over 
the wavelength range 0.95 to 1.1 /xm, which rough ly cor- 
responds to the Y band (Hillen brand et aL '2002"). The 
spectra of the L dwarfs arc dominated by absorption fea- 
tures arising from the 0—0 band of the F ^A — X ^A 
system of FeH. The strongest FeH feature is the band- 
head at 0.9896 /itm but weak FeH features are seen in 
the spectra of the L dwarfs throughout this wavelength 
range. The 0-0 band of the A *H - X '^E" system of 
VO centered at ~1.06 /im is also present in the spectra of 
the early-type L dwarfs and the FeH bandhead at 0.9896 
/xm is weakly present in the spectra of SDSS J1254— 0122 
(T2) and 2MASS 0559-1404 (T4.5). Finally, H2O and 
CH4 bands at A > 1.09 /im emerge at the L/T transition 
and strengthen through the T sequence. 

The synthetic spectra fit the overall shape of the Y- 
band spectra of the L and T dwarfs reasonably well. In 
particular, they match the strengths of the H2O and CH4 
bands in the T dwarf spectra at A > 1.08 /im. How- 
ever there are a number of mismatches worth noting. 
Firstly, the FeH bandhead in the synthetic spectra is 
systematically weaker than in the L dwarf data. Sec- 
ondly, all of the synthetic spectra exhibit the 0—1 band- 
head of the A 6e+ - X 6S+ system of CrH at 0.9969 
/im. This bandh ead is often identified i n the spectra 
of L dwar fs (e.g., iKirkpatrick et al.lll999t iKendall et al.l 
|2006[ ) but ICushing et all ()2003f ) questioned these iden- 
tifications given that many of the absorption features 
at this wavelength in _R=20 00 sp ectra can be associated 
with FeH. lReiners fc Basril (|2006[ ) recently found that no 
CrH absorption features can be identified in high res- 
olution (i?=31,000) spectra of M and L dwarfs which 
indicates that CrH is not a significant absorber at these 
wavelengths. Therefore the presence of this bandhead 
in the synthetic spectra is inconsistent with the obser- 
vations. The 0—1 ba ndhead of the (p system of TiO 
(jGalehouse et al.lfl980( l at ~1.104 /im is also too strong 
in all of the model spectra. The T dwarf models ex- 
hibit a broad, weak VO band centered at ~1.06 /xm. 
Observationally, this band peaks in strength at a spec- 
tral type of ^L5 and is absent by a spectral type of L8 
(jCushing et al.|[2005l ) so the presence of this band in the 
models is also inconsistent with the observations. Finally, 
the best fitting model spectrum underestimates the peak 
flux of the T4.5 dwarf 2MASS 0559-1404 because of the 
presence of the condensate opacity (/scd=4). A cloud- 
less model would match the overall shape of the Y band 
better (see J jiXTjl . 

4.3.2. The J Band 

Figure [5] shows the spectra and best fitting model spec- 
tra from 1.1 to 1.34 /im. This wavelength range contains 



the most prominent atomic features in the infrared spec- 
tra of L and T dwarfs, namely the two K I doublets at 
^1.17 /im and ^1.24 /xm, and the Na I doublet at 1.14 
/im. The 0-1 and 1-2 bands of the F '^A - X *A sys- 
tem of FeH are also present in the J-band spectra of the 
L dwarfs. These two bands exhibit bandheads at 1.1939 
and 1.2389 /im, respectively, along with numerous addi- 
tional weak features from 1.2 to 1.3 /im. The Sv^ and 
i/2 + 2i>'3 bands of CH4 centered at ^1.15 and '-^1.4 /im 
are also seen in the spectra of the T dwarfs. 

The model spectra reproduce the depths of the H2O 
bands at A >1.32 /im, the onset of the CH4 absorption in 
the spectra of the T dwarfs, and the weakening and even- 
tual disappearance of the Na I doublet at 1.14 /im with 
later spectral type. However the model spectra do a poor 
job of matching the peak flux levels at ^1.3 /im of 2MASS 
1439-M929 (LI), 2MASS 1507-1627 (L5), and 2MASS 
0559-1404 (T4.5). In addition, the bandheads and other 
weak FeH features observed in the J-band spectra of the 
L dwarfs are completely absent from the models. This 
behavior mirrors that in the K-band for which the FeH 
band is systematically weaker in the models than in the 
data. The cross-sections of the 0—1 and 1—2 /im FeH 
bands are an order of magnitude smal ler than those of 
the 0-0 FeH band (|DuHck et all 120031) so it is not sur- 
prising that they are either very weak or absent in the 
synthetic spectra. 

The models also are a poor match to the early-type 
L dwarfs near the 1.177 /im K I line since the red wing 
appears much too broad given the observations. This 
apparent line broadening in the model spectra is actu- 
ally due to th e presence of the 0— 1 band of the A— X 
system of VO (|Cheung et al.lll982l ). Although this band 
has been identi fied in the spectra of y oung, low-gravity 
brown dwarfs (jMcGovern et al.ll2004f) . it has not been 
found in the spectra of field L and T dwarfs and therefore 
its presence in the model spectrum implies the opacity 
of this VO band may be in error. 

4.3.3. The H Band 

Figure [9] shows the spectra and best fitting synthetic 
spectra over the 1.4 to 1.8 /im wavelength range. The L 
dwarf spectra are shaped primarily by H2O absorption 
but also exhibit numerous weak absorption features, in- 
cluding three bandheads, from 1.59 to 1.74 /im that arise 
from the E '^IV - A -^Tl system of FeH. The 2v:i and 
2i^2 + ^3 bands of CH4 appear at the L/T transition and 
strengthen through the T sequence. As explained in ij4.1i 
we do not include the 1.58 to 1.75 /im wavelength range 
in the fits because the FeH band lacks a line list and the 
CH4 line list is incomplete. 

The model spectra fit the overall shape of the data 
reasonably well at these wavelengths. In particular, the 
detailed structure of the H2O bands at the short and long 
wavelength ends of the H band are well matched by the 
model spectra. The best agreement between the data 
and model spectra at these wavelengths are obtained for 
DENIS 0255-4700 (L8) and SDSS 1254-0122 (T2) be- 
cause the FeH band is either weak or absent from the 
data and only weak CH4 absorption is present. This in- 
dicates that the H2O opacity tables are quite good in the 
H band. 

4.3.4. The K Band 
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Figure [10] shows the data and best fitting synthetic 
spectra over the 1.9 to 2.4 /im wavelength range. The 
spectra of L and T dwarfs are shaped primarily by H2O, 
CO, and CH4 bands at these wavelengths. At the earliest 
spectral types, CO is the dominant carbon bearing gas 
and therefore the Az^ = +2 CO overtone bands at A > 
2.29 /Ltm are strong. These CO bands weaken and disap- 
pear with later spectral type as CH4 becomes the dom- 
inant carbon-bearing gas. Additionally there is a Na I 
doublet at 2.26 /im present in the early- type L dwarfs. 

With the exception of 2MASS 1507-1627 (L5) and 
2MASS 0559-1404 (T4.5), the models are a poor fit 
to the i^-band spectra. The models underestimate 
the if -band flux levels of the 2MASS J1439-hl929 (LI) 
and 2MASS 0036-1-1821 (L3.5) and overestimate the 
strengths of the CH4 band in DENIS 0255-4700 (L8) 
and SDSS J1254-0122 (T2). In addition, a weak CH4 
feature can be seen in the synthetic spectrum of 2MASS 
0036-^1821 (L3.5) at '-^2.2 /itm that is not observed in 
spectra of ultracool dwarfs until a spectral type of ~L8. 
The presence of this CII4 feature in the model spec- 
trum of 2MASS 0036-1-1821 is a result of the fact that 
/scd=3 for this model. For a comparison, the best fitting 
model for 2MASS 1507-1627 (L5) has the same Teff, but 
/sod=2, and shows no CH4 feature (the change in g makes 
little difference in the strength of the CH4 band at this 
Tes). The /sed=3 model is cooler in the CH4 band for- 
mation region resulting in a higher CH4 abundance and 
a deeper absorption band. 

As described in iJ31 H2 CIA is an important pressure- 
sensitive opacity source in the atmospheres of ultracool 
dwarfs. Compared to the absorption of other molecules, 
(primarily H2O and CH4), H2 CIA is most significant in 
the K band. A lower matallicity decreases the opacity of 
H2O and CII4 relative to the II2 CIA continuum, which 
results in stronger emission in the Y , J and H bands with 
respect to the K band (JB urgasse r et al.ll2002al) . In previ- 
ous work (jSaumon et al.,.2006, .20071 ). we found that our 
solar metallicity models systematically underestimated 
the iiT-band fluxes of late-type T dwarfs while models 
with [Fe/H] = -f 0.3 fit the K-hand peak weU. Therefore 
the poor match between the data and model spectra in 
the K band may indicate that the metallicities of the L 
and T dwarfs in our sample differ somewhat form the 
solar value employed in our chemical equilibrium calcu- 
lations. Another possibility that explains the mismatch 
between the spectra DENIS 0255-4700 (L8) and SDSS 
1254-0122 (T2) and the best fitting models is the cloud 
properties since the cloud sedimentation efficiency fsed 
changes from 2 to 3 to 4 from L8 to T4.5. 

4.3.5. The L Band 

Figure [TT] shows the spectra and best fitting synthetic 
spectra over the 3.0 to 4.1 ^m wavelength range. The 
spectra are dominated by weak H2O features and the v^ 
fundamental band of CH4. The Q-branch (~3.3 /im) is 
the first CII4 feature to emerge in the spectra of the L 
dwarfs and does so at a spectral type of '^L5. With in- 
creasing spectral type, the P and the R branches emerge 
and strengthen through the spectral sequence until the 
band is saturated in the late- type T dwarfs. 

The model spectra do not match the data particularly 
well in the L band. The most egregious mismatch is 
2MASS 0036-^1821 (L3.5) for which the model shows a 



very strong CH4 band that is absent in the data. The 
presence of the CH4 band is due to the fact that the best 
fitting model has /sod=3, as in the K band (see H4.3.5p . 
The average flux density levels of the models also ap- 
pear low for 2MASS 1507-1627 (L5), DENIS 0255-4700 
(L8), SDSS J1254-0122 (T2), and 2MASS 0559-1404 
(T4.5). In addition, the model for 2MASS J0559-1404 
does not match the peak of the observed spectrum at 
~4.1 /im. Generally, in objects where the CH4 absorp- 
tion is too strong in the L band, a corresponding mis- 
match occurs in the K and IRS/SL bands (L8 and T2). 
The line lists for both H2O and GH4 are reasonably com- 
plete at these wavelengths so the most likely explanation 
for the mismatch is that the region of the model atmo- 
sphere from which the K and L band and IRS/SL flux 
emerges is too cool, possibly as a result of the condensate 
cloud model. 

4.3.6. IRS/SL 

Figure [T^ shows the data and best fitting models over 
the 5.5 to 14.5 fim waveleng th range covered by the 
Short-Low module of the IRS (iHouck et al.l[200l . H2O 
absorption features arising from the 1^2 fundamental band 
and 2i/2 ~ 1^2 overtone band are present throughout this 
wavelength range; the break in the spectra at 6.5 /im is 
an H2O feature. The v^ fundamental band of CH4 cen- 
tered at ~7.65 fim emerges in the spectra of late-type L 
dwarfs and the 1^2 fundamental band of NH3 centered at 
^^10.5 /im appears in the spectra of the T dwarfs. 

With the exception of 2MASS 0036-^1821 (L3.5), the 
synthetic spectra match the data reasonably well. In 
particular, the strengths of the H2O band at 6.5 /xm, 
and the CH4 band at 7.8 /im is well matched by the 
model spectra. The flux level of the synthetic spectrum 
of DENIS 0255-4700 (L8) is low however, and there 
is a slight mismatch between the model spectrum and 
2MASS 1507- 1627 (L5) at ~9 /im that we discuss fur- 
ther in ij4.4.5l The mismatch between the model and the 
spectrum of 2MASS 0036- H821 is d ue to t he fact that 
the model has /sod =3 (see ^4.3.41 and ^4X5)) . 

4.4. Fits to Individual Photometric Bands 

The L and T dwarf spectra studied herein cover a 
broad range of wavelengths while most published spec- 
tra cover only the red o ptical a nd near-i nfrared wave- 
lengths (e.g., ICruz et al . 200^ iKnapp et al. 20Q3) or 
even a single photometric band (e .g. 
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120051 iMetchev fc Hillenbrand! I2006D . It is therefore ap- 
propriate to investigate how well the atmospheric param- 
eters of L and T dwarfs can be determined by fitting L 
and T dwarf spectra over narrower wavelengths ranges. 
In addition, we can also estimate the systematic errors in 
both the models and spectral fitting procedure by com- 
paring how well the models fit the spectral features over 
narrower wavelength ranges as compared to the global 
0.95-14.5 /im fits. 

We therefore fit the data over the wavelength ranges 
described in §4.31 as well as over the 0.95—2.5 fim wave- 
length range. For the IRS/SL and Y-, J-, K-, and L- 
band fits, we set Wi — 1 for all i. Table [5] lists the results 
of the fits and Figure [T21 which is similar to the lower 
panel of Figure shows the corresponding T^g values for 
each dwarf. Figures [Tll-fTBl and [THH^D] show the spectra 
and best fitting models {blue) along with the best fitting 
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models obtained by fitting the 0.95—14.5 /im {red, same 
as Fig. [6]). 

As expected, the model spectra fit the data much bet- 
ter over these narrower wavelength ranges. In most cases, 
remarkably good fits are obtained. However, the T^e 
values derived from narrower wavelength ranges do not 
agree with those derived by fitting the 0.95—14.5 fim 
SEDs. Particularly poor Tcs values are derived by fit- 
ting the L band and IRS/SL spectra of the early- to 
mid- type L dwarfs (see Figure [T5|) . Below, we describe 
how the model spectra fit the data over the narrower 
wavelength ranges in more detail. Since the majority of 
the H band is not included in the fitting process, (see 
M.l\i . we do not present a detailed description of that 
wavelength region below. 

4.4.1. Y Band Fits 

Figure [14] shows the spectra and best fitting model 
spectra over the Y band. The fits reproduce the over- 
all shape of the Y band data considerably better than 
0.95—14.5 /im fits for all spectral types. The remaining 
descrepencies highlight problems with specific molecu- 
lar bands. The model FeH bandheads at 0.9896 /im are 
still too weak, the 0.8611 /im CrH bandheads are still 
present, and the TiO bandheads at ~ 1.104 /im are still 
too strong. Indeed the CrH bandhead is present in all the 
model spectra shown in Figure [Mj The strength of the 
~1.06 /im VO band is, however, much better matched 
when fitting the Y band data alone. The peak flux of 
2MASS 0559-1404 (T4.5) is also better matched when 
fitting the Y band data alone but the synthetic spectrum 
does exhibit the VO band that is absent from the data. 

As shown in Figure fT3| six of the ni ne dwarfs have de- 
rived Tpff values that fall within the iGolimowski et all 
()2004( ) effective temperature range. However the derived 
Teff values for the early- to mid- type L dwarfs (LI to L5) 
differ by no more than 200 K while the IGolimowski et al.l 
()2004[) 3-Gyr effective temperatures and the TgB values 
derived by fitting the 0.95—14.5 /im spectra vary by 500 
K over the same spectral type range which indicates the 
Y band is not a sensitive indicator of effective tempera- 
ture at these spectral types. 

4.4.2. J Band Fits 

Figure [15] shows the spectra and best fitting models 
over the 1.1 to 1.34 /im wavelength range. Fits to the J- 
band spectra alone significantly improve for those dwarfs 
where the global fit was rather poor, i.e., for 2MASS 
1439+1929 (LI), 2MASS 1507-1627 (L5), and 2MASS 
0559-1404 (T4.5). The FeH features are stiU absent from 
the models in the early- to mid-type L dwarfs, but over- 
all, very good fits in the J band can be obtained for the 
whole spectral sequence. Except for 2MASS 1439-f 1929 
(LI) and 2MASS 0036-^1821 (L3.5), the T^s values ob- 
tained with the J band arc consistent with those derived 
from fitting the fuU SED. 

4.4.3. K Band Fits 

Figure [16] shows the spectra and best fitting models 
over the 1.9 to 2.4 /im wavelength range. The improve- 
ment in the quality of the fits is substantial when the 
-fC-band spectra are fitted alone. Indeed the agreement 
between the synthetic spectra and data is at its high- 
est at these wavelengths; the detailed structure of the 



H2O, CO, and CH4 absorption features throughout this 
wavelength range are almost perfectly reproduced by the 
models (see Figure [T71) . However the improvement in the 
fits comes from a change in T^e of up to 300 K and t he K - 
band Teff values show more scatter. Bla ke et al.l ()2007t ) 
have fitted high resolution {R « 50000) 2.296-2.308 /xm 
spectra of nine L dwarfs spanning spectral types LO to L6 
using our model atmospheres and also find good agree- 
ment between the data and models. 

4.4.4. L Band Fits 

Figure [TH] shows the spectra and best fitting mod- 
els over the 3.0 to 4.1 /im wavelength range. There is 
a marked improvement in the quality of the fits when 
the L-band spectra are fitted separately. The early- to 
mid-type L dwarfs no longer show the deep Q-branch 
CH4 absorption and the CH4 band in DENIS 0255-4700 
(L8), SDSS J1254-0122 (T2), and 2MASS 0559-1404 
(T4.5) is well matched by the synthetic spectra. The 
model spectrum now fits the peak flux in the spectrum 
of 2MASS 0559-1404 at 4.1 /im well. The derived T^s 
values are systematically hotter than those we obtained 
from the other bands by up to 500 K for all the dwarfs 
except 2MASS 0559-1404 (Fig. [TS]) which indicates that 
the L band is a poor spectral region from which to de- 
rive atmospheric parameters of L and T dwarfs. The low 
value of /mc for some of the dwarfs in our sample (Table 
[5]) reinforces this conclusion. 

4.4.5. IRS/SL Fits 

Figure [19] shows the data and best fitting models over 
the 5.5 to 14.5 /zm wavelength range. There is also a 
marked improvement in the fits over this wavelength 
range since the strengths of the H2O, CH4, and NH3 
bands are we ll reproduced by the model spectra. As 
first noted bv lCushing et al.] (|2006D . there is a mismatch 
between the data and the models from 9 to 11 /xm in 
the spectra of the mid-type L dwarfs. This mismatch is 
most prominen t in the spectrum of 2MASS 2224-0148 
(c.f.. Fig. 9, ICushin g et al.l [20061 ) but can also been 
seen weakly in 2 MASS 0036-H 821 (L3.5) and 2MASS 
1507—1627 (L5). ICushing et al.l have tentatively identi- 
fied this feature as arising from the Si-0 stretching mode 
of small silicate g rains that give rise to the classic 10 /im 
silicate feature. iHelling et al.l ()2006f ) suggest that the 
absorption feature could arise from quartz (Si02) grains, 
although this is not expected from equilibrium chemistry. 

The derived Tcs values for the early- to mid-type L 
dwarfs (Fig. [0)1 are well above the globally determined 
values (Fig. [5]) by up to 700 K. For four dwarfs, the best 
fit is found in the high- Toff high-g corner of our model 
grid (Teff =2400 K, log(/=5.5) and should be considered 
with caution. The IRS/SL spectral range gives Teff val- 
ues that behave much like those obtained by fitting the L 
band. Both the L and IRS/SL bands probe regions of the 
atmosphere with the the lowest pressures and tempera- 
tures. The tendency of the L and IRS/SL band fits to 
produce high Teff values combined with the fact that the 
global 0.9—14.5 /im fits tend to over predict the depth of 
the 3.3 and 7.8 /im methane bands implies that the tops 
of the model atmospheres may be too cold. This could 
be simply due to flux redistribution caused by inadequate 
CH4 opacities at A < 1.7 /im or a different vertical distri- 
bution of opacity within the cloud. On the other hand. 
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atmosphere models also tend to under-predict the tem- 
peratures of the middle a nd upper atmospheres o f solar 
system giant planets (^e.g.. lMarlev &: McKavlfTOOgf ) . even 
when accounting for the effects of incident radiation and 
photochemistry. One possible explanation for this dis- 
crepancy is that upwardly propagating waves launched 
within the convection region break and deposit their en- 
ergy in the lower pressure regions of the atmosphere (e.g., 
the discussion in lYoung et al.l (|2005f ) for Jupiter). An ex- 
ploration of the possible effects of such waves in ultracool 
dwar f atmospher e s wou ld be o f interest. 

Sa umon et all (|2006L 120071 ) have found that the 
IRS/SL spectra of Gl 570D (T7.5) and 2MASS 
J04151954-0935066 (T7) can be modeled adequately 
only by reducing the abundance of NH3 by about an 
order magnitude from that predicted by chemical equi- 
librium. This is attributed the vertical motions of the gas 
within the atmospheres which prevents the abundances 
of N2 a nd NH3 from reaching th eir expected equilibrium 
values (iLodders fc Feglev '2002*) . The models presented 
herein do not include the effects of non-equilibrium chem- 
istry. We note that an equilibrium model fits the IRS/SL 
spectrum of 2MASS 0559-1404 (T4.5) well. This is a 
consequence of its higher T^s, which decreases the effect 
of vertical transport on the chemistry of NH3. Neverthe- 
less, a non-equilibrium model provides a better fit to the 
NH3 features of the IRS Sh ort High spectrum ( RfvGOO) 
of this object at a 2-a level (jMainzer et al.ll2007f ). 

4.4.6. 0.95-2.5 ^im Fits 

Figure [20] shows the data and best fitting models over 
the 0.95 to 2.5 ^m wavelength range. With the excep- 
tions of 2MASS 1439-H1929 (LI) and SDSS 1254-0122 
(T2), fitting the near-infrared wavelengths alone does 
not dramatically improve the quality of the fits. In- 
deed the same model ([1700/5.5/3]) fits both the near- 
infrared and 0.95-14.5 ^lm spectra of 2MASS 0036+1821 
(L3.5). StiU, the depths of the H2O and GH4 bands 
in the remaining dwarfs are better matched when fit- 
ting the near-infrared spectra alone. Note also that 
the near-infrared spectra of 2MASS 0036+1821 (L3.5) 
and 2MASS 1507-1627 (L5) are fit by the same model 
(1700/5. 5/f3). The derived T^s values do not match the 
0.95—14.5 /im T^e values well, particularly for the early- 
to mid-type L dwarfs. This implies that spectra at A > 
2.5 fim are necessary in order to derive Teg values for 
early- to mid-type L dwarfs. This is unfortunate since 
most published spectra of L dwarfs rarely extend beyond 
2.5 /im. 

4.5. The L/T Transition 

The transition from the L to T dwarfs is marked by 
a rapid change in near-infrared color from the very red 
late-type L dwarfs (J — K^2.b) to the very blue T 
dwarfs (J — K'^—1) an d an inc rease in the emergent 
flux in the 1 txai r egion (jDahn et al. 2003; iTinnev et al.l 
120031 IVrbaet al.l l2004f). The latter trend was unex- 
pected because theoretical Lboi declines as a brown dwarf 
cools. The 0.9896 /im FeH bandhead also strengthens 
in the spectra of the early- to mid-type T dwarfs after 
weakening considerably in the spectra of the late-typ e 
L dwarfs fBurgass er et al.ll2002bl : iMcLean et al.ll2003f ). 
Furthermore, the transition occurs ove r a very small Teg 
range, perhaps only 100 to 200 K (jGolimowski et ahl 



|2004[). All current atmospheric and evolutionary mod- 
els (jMarlev et al.ll2002t lTsuiill2005t iBurrows et al]l2006f l 
fail to reproduce these observations. 

A number of possible solutions have been proposed 
including the break up of t he condensate cloud d ecks 
(Ackerman fc MarlevI [2OOII : iBurgasser et al.l l2002bf ). a 
rapid c hange in the efficie ncy of condensate sedimen- 
tation (iKnapp et"all I2004D . or even "crypto-binarity" 
(jBurrows et al.l[2006[ ) . Unresolved binaries do play a role 
in exaggerating the magnit ude of the 1 /xm b rightening 
in the field population (Burgasser et al .ll2006bl: Burgassen 
l2007t iLiu et al1l2006( l but the brightening appears to be 
intrinsic to the evolution of brown dwarf atmospheres 
since it has also been observed in binary systems com- 
posed of early-type and mid-type T dwarfs. In these 
systems, whose components presumably have the same 
metallicity and are coeval, the compone nt with the later 
spectral type is brigh ter in the J band (jBurgasser et al.l 
[20M[LhreraII[200m . 

Although the models involving changes in the conden- 
sate cloud properties are both consistent with the color 
magnitude diagrams of field L and T dwarfs at the L/T 
transition (jBurgasser et all l2002bl : iKnapp et aLll2004( l. 
the model spectra have yet to be directly compared to 
spectroscopic observations. Indeed an error in the TiO 
line list used in a previous generation of our models (see 
SjS]) produced emergent spectra with TiO bands at near- 
infrared wavelengths that were much too strong, but 
whose J- and K-hand magnitudes nevertheless matched 
the colors of L dwarfs relatively well. Unfortunately, self- 
consistent models that include the break-up of the con- 
densate clouds decks are not yet available so we can only 
compare models with changes in the sedimentation effi- 
ciency /scd- 

Four of the dwarfs in our sample span the L/T tran- 
sition: 2MASS 0825+2115 (L7.5), DENIS 0255-4700 
(L8), SDSS J1254-0122 (T2), and 2MASS 0559-1404 
(T4.5). As can be seen in Figure [S] and Table [31 the 
effective temperatures derived from fitting the 0.9—14.5 
/im spectra of these four dwarfs decrease slightly from 
1400 K to 1200 K and the clouds become gradually thin- 
ner (/scd increases) toward later spectral types. Our fits 
of the SEDs of late-type L and early-type T dwarfs there- 
fore support the concept of a rapid decrease in the overall 
cl oud opacity over a narro w range in Toff. 

IBurgasser et al.l (|2002bD argue that the resurgence of 
the 0.9896 /im FeH bandhead i n the spectra of ear ly- 
to mid-type T dwarfs (see also iMcLean et al.l l2003l ) is 
consistent with the cloud breakup hypothesis. In the ab- 
sence of such a breakup, radiation at ~1 /im only emerges 
from atmospheric layers above the Fe cloud deck in which 
gaseous FeH has almost completely been depleted into 
solid or liquid Fe. If holes are present in the cloud decks, 
radiation from atmospheric layers below the Fe cloud 
deck where FeH has not been depleted can escape and as 
a result, the FeH bandhead is observed in the spectra. 

Figure [21] shows a sequence of models with rcff=1400 
K, log 5=5.0, and f^^^—l^ 2, 3, 4, and nc. The strength of 
the FeH bandhead also increases in strength with increas- 
ing sedimentation efficiency. For a fixed P /T profile, the 
FeH column density above the Fe cloud base (T ~ 2000 
K) is set by the vapor pressure of Fe and thus is indepen- 
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dent of both nietallicity and /scd^^- In the /sod=l niodel, 
the T - 2/3 level is high in the cloud (T w 1500 K) but 
as /sod increases, the Fe cloud becomes thinner and the 
r ^ 2/3 level is reached deeper in the atmosphere. As 
a result, the column density of FeH that is observed in- 
creases resulting in an increase in the band strength of 
FeH. Therefore a change in sedimentation efhciency at 
a fixed Tcs can also explain the resurgence of the FeH 
bandhead across the L/T transition. 



4.6. Very Red L Dwarfs 

As can be seen in Figure (TJ there is significant scat- 
ter (^0.5—1 mag) in the near-infrared colors of L and 
T dwarfs at a given spectral type. It has been sug- 
gested that this spread ari ses from variations in the 
condensate c loud properties (iKnapp et al.ll2004 ) or sur- 
face gravity (jBurrows et al.l 120061 ), altho ugh unresolved 
binar ies also contribute to the scatter (jLiu fc Leggettl 
l2005f) . In our models, both surface gravity and cloud 
sedimentation efficiency affect the near-infrared col- 
ors, but the latter dominates for the mid- to late- 
type L dwarfs. Two pairs of objects in our sam- 
ple, 2MASS 2224-1058/2MASS 1507-1627 and 2MASS 
J0825-h2115/DENIS J0255-4700, have identical spectral 
types within the errors but markedly different J — Kg 
colors: 2.05/1.52 and 2.07/1.69, respectively. As a re- 
sult, they are ideal for studying the dependence of sur- 
face gravity and the condensate cloud properties on the 
colors of L dwarfs at a given spectral type. 

The best fitting model spectra over the 0.95 to 14.5 
/im wavelength range for 2MASS 2224-0158 and 2MASS 
1507-1627 have [1700/4.5/1] and [1700/4.5/2], respec- 
tively. The two objects have the same log 5 and T^s 
values, but the redder object, 2MASS 2224-0158 has 
an /sod value of 1 while the bluer 2MASS J1507-1627 
has /sod=2. Condensate clouds with /sod=l are thicker 
and have smaller particles and therefore contribute more 
opacity to the atmos phere (which res u lt in redder J — 
Ks colors (see !j3])). iBurgasser et al.l (|2007[ ) have also 
noted that the broad absorption feature s een in the IRS 
spectrum of 2MASS J2224-0158, which 'Cushi ng etlHI 
(|2006) have tentatively identified as arising from small 
silicate grains in the condensate clouds (see iJ4.4.5p . is 
much weaker in the IRS spectrum of the bluer 2MASS 
1507—1627. This is consistent with the condensate 
clouds being responsible for the very red color of 2MASS 
2224-0158. 

Similarly, the best fitting model spectra over the 0.95 
to 14.5 urn. wavelength range for 2MASS 0825-^2115 and 
DENIS 0255-4700 have [1400/4.5/1] and [1400/5.5/2], 
respectively. These two objects have identical T^s and 
different logg values but as above, the redder object, 
2MASS 0825-1-2115, has /sed=l- These latter two objects 
show how both surface gravity and the condensate clouds 
play a role in controlling the colors of L dwarfs. 

^"^ The formation of FeH gas can be written as Fe (gas) -|- 0.5 H2 
= FeH (gas). If Fe is condensed, the partial pressure of Fe (gas) 
is fixed by the vapor pressure of the iron and is independent of 
the total amount of condensed phase. Since the partial pressure of 
H2 is essentially a constant (at a given total pressure), the partial 
pressure of FeH is also fixed. Therefore once Fe condenses, the 
partial pressure, and thus column density, of FeH above the Fe 
cloud base is independent of both metallicity and /sod for ^ fixed 
P/T profile. 



5. SYSTEMATIC UNCERTAINTIES 

Both the G statistic and the Monte Carlo simulations 
used in this analysis explicitly account for the random 
errors in the data. However it is important to also un- 
derstand any systematic errors that can affect the anal- 
ysis. Below, we discuss two such systematic errors, the 
absolute flux calibration process and the choice of the 
weights in the G statistic. 

5.1. Absolute Flux Calibration 

The near-infrared spectra used in this work were ob- 
tained with SpeX on the NASA IRTF. The 0.8-2.5 yum 
wavelength range is covered in six cross-dispersed orders, 
and as a result the spectra from the various orders are 
merged into a single spectrum by scaling the spectra to 
ensure that the flux density levels match in regions where 
the spectra overlap in wavelength. Unfortunately there is 
no wavelength overlap between the H and K bands and 
therefore no scaling could be performed. The merged 
0.8—2.5 /xm spectra were then absolutely flux calibrated 
by deriving a single scale factor based o n the 2MASS 
JHKs photometry. ICushing et al.l ()2005[ ) showed that, 
on average, the synthetic colors {J — H, H — Kg, and 
J — Ks) of their entire sample of M, L, and T dwarfs 
match the observed 2MASS colors within the errors. 
However in some cases, we find deviations between the 
synthetic and observed colors of up to 2 a. 

In order to quantify any effect that variations in the 
absolute flux calibration of the near-infrared spectra have 
on our analysis, we have absolutely flux calibrated each 
of the SpeX cross-dispersed spectra separately using both 
2MASS photometry a nd Mauna Kea Observ atories Near- 
Infrared (MKO-NIR ; iTokunaga erani2002f) photometry 
(jKnapp et al.ll20M )^^. Two of the dwarfs in our sample, 
2MASS 1506-^1321 and DENIS J0255-4700, have not 
been observed in the MKO-NIR system to date. The 
near-infrared, L band, and IRS/SL spectra are then re- 
combined into a single 0.95—14.5 /.tm spectrum and the 
entire analysis is performed again. 

Figure [H] summarizes the differences between the at- 
mospheric parameters derived by fitting the spectra with 
the different absolute flux calibration schemes over the 
0.95—14.5 /im wavelength range. Flux calibrating the 
spectra with 2MASS photometry leaves the derived pa- 
rameters almost unchanged while the derived parame- 
ters for three objects flux calibrated with the MKO-NIR 
photometry change by 100 K in T^s and 0.5 dex in logg. 
Overall, there is little change in the derived values of 
Tcff', log 5, and fscd which indicates the absolute flux cal- 
ibration does not severely affect the parameters derived 
by fitting the 0.95-14.5 /xm spectra (Figure!^. The 
parameters obtained by fitting the 0.95—2.5 /im spectra 
alone are more sensitive to variations in the absolute flux 
calibration. Differences of 0.5—1 dex in logg, 100 K in 
Tcff, and 1 in /sod in the best fitting model spectra are 
apparent. The range of possible models (i.e., the error 
bars) also grows because the SpeX cross-dispersed orders 
(J, iJ, K) are absolutely flux calibrated separately and 
thus are varied separately in the Monte Carlo simulation. 

^^ The 2MASS J band filter spans two spectral orders and there- 
fore we have scaled these two orders to the same flux density levels 
by eye. The two orders are then absolutely flux calibrated together 
using J-band photometry. 
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5.2. Weights 

As discussed in i i4.11 we weight each pixel in the spec- 
trum by its width (wi = AA^) to remove a bias introduced 
into the G statistic due to the large variation in wave- 
length sampling of the data. Our choice of weights, how- 
ever, is completely arbitrary. In order to estimate any 
systematic error in the derived parameters introduced 
by our choice of weighting scheme, we have devised a 
second weighting scheme as follows. 

We divide the 0.95—14.5 /im wavelength range into 
bins equally-spaced in In A. Each pixel that falls in bin 
k receives a weight inversely proportional to the total 
number of pixels contained in bin k. For example, if bin 
k contains four pixels, then 'Wi=0.25 for each pixel in bin 
k. The Spitzer data has the lowest resolving power of 
all the data (i? = A/AA sa 100) so we used a bin size of 
6\nX = 0.01 which corresponds to i?=100. This alterna- 
tive weighting scheme effectively mimics smoothing the 
spectra to i?=100 and resampling them such that each 
pixel corresponds to a single resolution element. 

Figure [21] summarizes the differences between the at- 
mospheric parameters derived by fitting the data using 
the two different weight schemes over the 0.95—14.5 fim 
wavelength range. The Toff values derived for the early- 
type L dwarfs are lower by 100 K and the derived logg 
values for two of the dwarfs change by 1 and 0.5. Based 
on these results, the choice of weighting scheme does not 
dramatically alter the results presented herein. 

6. CONCLUSIONS 

We have performed the first fits of the 0.95— 14.5 /im 
SEDs of nine field ultracool dwarfs spanning spectral 
types LI to T4.5 to obtain effective temperatures, grav- 
ities and values for the cloud sedimentation parameter. 
The wide wavelength coverage of the data allows us to 
easily discriminate between synthetic spectra with differ- 
ent atmospheric parameters because most physical and 
chemical processes that occur in ultracool dwarf atmo- 
spheres leave a signature in one or more of the wavelength 
regions covered by the data. 

In our study of field ultracool dwarfs we have focused 
on the most significant physical parameters, ignoring the 
more subtle effects of metallicity variations from the so- 
lar value as well as departures from chemical equilib- 
rium caused b y verti cal mixing in the atmosphere (e.g. 
ISaumon et al.l (j2006f )). The spectroscopic data come 
from a variety of instruments and we have estimated the 
uncertainties in the derived parameters that arise from 
the absolute fiux calibrations of the various spectra that 
form the SEDs. Although some wavelength ranges are 
poorly fitted for some objects, undoubtedly due to re- 
maining inadequacies in the models and, to a smaller 
extent, to possible variations in metallicity amongst the 
dwarfs in our sample, overall the model spectra fit the 
data well. 

Given the natural scatter in the observational proper- 
ties among the objects in our sample and the uncertain- 
ties inherent in our procedure, we can identify only gen- 
eral trends in the fitted parameters as a function of spec- 
tral type. The effective temperatures decrease steadily 
through the LI to T4.5 spectral type s and agree with the 
determinations of iGolimowski et al.l (|2004). Our results 
also confirm that T^g is nearly constant at the L/T tran- 
sition, decreasing by only ^ 200 K from spectral types 



L7.5 to T4.5. 

Two topics of great interest are the evolution of the 
cloud properties as a function of spectral type and the 
distribution of cloud properties within each spectral sub- 
class. Our work provides a first, although somewhat 
blurry glimpse of this important aspect of ultracool dwarf 
atmospheres. We find that the cloud properties, as en- 
capsulated in the most imp ortant free parameter of the 
lAckerman fc MarlevI (|2001h condensation cloud model, 
/scd, vary significantly among L dwarfs, ranging from 
very thick clouds (/sod = 1) to relatively thin clouds 
(/sod — 3) with no particular trend with spectral type. 
The two L dwarfs in our sample with very red near- 
infrared colors have best fitting model spectra with thick 
clouds. Finally the fits to the two T dwarfs in our sample 
(T2 and T4.5) suggest that the clouds becomes thinner 
in this spectral class, in agreement with previous studies. 

Fitting individual spectral bands, we almost always ob- 
tain excellent fits to the data, in some cases down to the 
fine structure of the H2O and GII4 bands. The Toff ob- 
tained in this fashion (Fig. I13p can show quite a bit of 
scatter compared to the values derived by fitting the full 
SED, typically by ^ 200 K and in the worst cases, up 
to 700 K. Nevertheless, the trend of decreasing Toff with 
spectral type is preserved in all bands, with few excep- 
tions. However, this exercise shows that in general. Toff 
determinations based on the analysis of narrow spectral 
ranges are not reliable for ultracool dwarfs. 

The data point to several directions for improving the 
models. First, improved input molecular opacities, par- 
ticularly for GrH, VO, FeH, and GII4 and the reson ance 
doublet of K I (7665, 7699 A) (lAllard et al.l[200l are 
required to improve fits to the early L and all of the T 
dwarfs, respectively. The second important area ripe for 
improvement remains the cloud model, in particular the 
particles' size distribution, the cloud's vertical profile, 
and the optical properties of the condensates. 

The accurate determination of the physical parame- 
ters of isolated field ultracool dwarfs therefore remains 
an elusive goal. Our method of SED fitting should be 
applied to those ultracool dwarfs with well constrained 
Toff (±20 K) and logg (±0.05 dex) values (derived from 
other methods) in order to better assess its internal un- 
certainties. Current ly, such a sample is limited to a few 
late-t ype T dwarfs (jSaumon et al.ll2006l 120071 : iLiu et all 
l2007f ) that are companions to main sequence stars with 
known distances, metallicities and well-constrained ages. 
Such studies should be extended to include earlier spec- 
tral types in order to firmly associate the observable 
properties of the L and T spectral sequence to funda- 
mental astrophysical parameters. 
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Fig. 1. — 2MASS J — H {top), H — Ks {middle), and J — Ks {bottom) colors as a function of spectral type for all of the L and T dwarfs 
in the DwarfArchives web site as of 2007—06—07. The dwarfs in our sample are shown in red. 
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Fig. 6.— The 0.95-14.5 ^lra spectra of 2MASS 1439+1929 (LI), 2MASS 0036+1821 (L3.5), 2MASS 1507-1627 (L5), DENIS 0255-4700 
(L8), SDSS 1254-0122 (T2), and 2MASS 0559-1404 {T4.5) are shown in black. The data have been normaUzed to unity at 1.27 /xm and 
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Fig. 7. — Same as Figure [6] except the data cover the Y band and were normalized at 1.05 fira. 
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Fig. 8. — Same as Figure [6] except the data cover the J band and were normalized at 1.25 fira. 
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Fig. 9. — Same as Figure [6] except the data cover the H band and were normahzed at 1.6 /ira. 



Model Comparisons 



25 



CO 

-*— ' 

CO 

o 
O 

+ 



OJ 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



H2O 



H2O 



L8 
(1400/5.5/f2) 




1.9 2.0 2.1 2.2 2.3 2.4 1.9 2.0 2.1 2.2 2.3 2.4 
Wavelength (|jm) Wavelength {[im) 



Fig. 10. — Same as Figure |6] except the data cover the K band and were normahzed at 2.1 jivci. 
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Fig. 11. — Same as Figure |6] except the data cover the L band and were normalized at 3.8 /im. 
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Fig. 12. — Same as Figure |6] except the data cover the IRS/SL wavelength range and were normaUzed at 12 jim. 
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Fig. 14.— The Y-band spectra of 2MASS 1439+1929 (LI), 2MASS 0036+1821 (L3.5), 2MASS 1507-1627 (L5), DENIS 0255-4700 
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Fig. 15. — Same as Figure [141 except the data cover the J band and were normalized at 1.25 fira. 
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Fig. 16. — Same as Figure [141 except the data cover the K band and were normalized at 2.1 lira. 
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the relative flux levels between the data and models. 
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Fig. 18. — Same as Figure [141 except the data cover the L band and were normalized at 3.8 fim. 
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Fig. 19. — Same as Figure [141 except the data cover the IRS/SL wavelength range and were normahzed at 12 fim. 
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Fig. 20. — Same as Figure [141 except the data cover the 0.95—2.5 fira wavelength range and were normahzed at 1.27 /ira. 
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Fig. 21. — Synthetic spectra smoothed to i?=2000 with Tcff=1400 K, logg=5.0, and /scd=l, 2, 3, 4, and the nc model. The depth of 
the FeH bandhead at 0.9896 fim increases as f^cd increases in value. 
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Fig. 22. — Effects of different absolute flux calibrations on the derived atmospheric parameters over the 0.95—14.5 fim wavelength range. 
The panels show the difference between the values derived after scaling the near-infrared spectra by a single scale factor and the values 
derived after scaling the J, H, and K band spectra separately using 2MASS (left) and MKO-NIR (right) photometry. 



38 



Gushing et al. 



I 



E 
o 



+ 1 - 
0- 

-1 

+2 

+1 [ 



M-"" 



"T — ■ — I — ■ — I — ' — I — ■ — I — ■ — I — ' — r- 



2MASS. . 






_l 1 I 1 I 1 I 1 I 1 I 1 L_ 




-1 ■ 



-2 



-| — ■ — I — ■ — I — ' — I — ■ — I — ■ — I — ' — r- 



LI... 



I 



_l 1 I 1 I 1 I 1 I 1 I 1 L_ 



"T ■ 1 ' 1 ■ 1 ' 1 ' 1 ■ T" 



l..li 



+200 




_ +100 


( 


» 


= 

1- 

^ -100 


-4 i i 


-200 






_l 1 I 1 I 1 I 1 I 1 I 1 l_ 



-| ■ 1 ' 1 ■ 1 ' 1 ' 1 ■ T" 



_l 1 I 1 I 1 I 1 I 1 I 1 l_ 



LO L4 L8 T2 LO L4 L8 T2 
Spectral Type 



Fig. 23. — Same as Figure l22l except the data were only fit from 0.95 to 2.5 iim. 
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Fig. 24. — Effects of different weigliting schemes on the derived atmospheric parameters over the 0.95—14.5 /.tm wavelength range. The 
panels show the difference between the values derived by setting the weights to Wi = AA^ and tlic values derived by setting the weights to 
Wi = 1/n where n is the number of pixels in bins of width 51nA=0.01. 
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TABLE 1 
The Properties of the L and T Dwarfs in Our Sample 



Object 


Spectral Type'' 


(mag) 


J-H 

(mag) 


H-K^ 




J-Ks 

(mag) 


(mas) 


Rcf. 




Optical 


intrared 


(mag) 




2MASS J14392836+1929149 


LI 


LI 


12.759±0.019 


0.718±0.027 


0.495±0.029 


1.213±0.029 


69.6±0.5 


1 


2MASS J15065441+1321060 


L3 




13.365±0.023 


0.985±0.031 


0.639±0.028 


1.624±0.030 






2MASS J00361617+1821104 


L3.5 


L4±l 


12.466±0.027 


0.878±0.040 


0.530±0.036 


1.408±0.034 


114.2±0.8 


1 


2MASS J22244381-0158521 


L4.5 


L3.5 


14.073±0.027 


1.255±0.037 


0.796±0.035 


2.051±0.035 


87.02±0.89 


1,2 


2MASS J15074769-1627386 


L5 


L5.5 


12.830±0.027 


0.935±0.036 


0.583±0.035 


1.518±0.037 


136.4±0.6 


1 


2MASS J08251968+2115521 


L7.5 


L6 


15.100±0.034 


1.308±0.047 


0.764±0.041 


2.072±0.043 


94.22±0.99 


1.2 


DENIS 025503.3-470049.0 


L8 


L9 


13.246±0.027 


1.042±0.036 


0.646±0.034 


1.688±0.036 


201.37±3.89 


3 


SDSS J125453.90-012247.4 


T2 


T2 


14.891±0.035 


0.801±0.043 


0.253±0.060 


1.054±0.064 


73.96±1.59 


2.4 


2MASS J05591914-1404488 


T5 


T4.5 


13.802±0.024 


0.123±0.050 


0.102±0.068 


0.225±0.057 


96.73±0.96 


1,2 




Jl|2002|l, 
rfs are fi 
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TABLE 2 
Spectral Characteristics of Ultracool Dwarfs in Our Sample 



Object 






R= 


A/AA 








Rcf. 




(0.6 


-0.9 fJ.m) 


(0.9-2.5 Mm) 


(3.0- 


4.1 /^m) 


(b.b- 


-14.5 fim) 




2MASS J1439+1929 




890 


2000 




940 




90 


1,2,2,3 


2MASS J1506+1321 




890 


2000 




940 




90 


4,2,2,3 


2MASS J0036+1821 




890 


2000 




940 




90 


5,2,2,3 


2MASS J2224-0158 




890 


2000 




425 




90 


6,2,2,3 


2MASS J1507-1627 




890 


2000 




940 




90 


5,2,2,3 


2MASS J0825+2115 




890 


1200 




210 




90 


6,2,2,3 


DENIS 0255-4700 




890 


2000 




210 




90 


7,2,2,3 


SDSS J1254-0122 




1200 


1200 




210 




90 


8,2,2,3 


2MASS J0559-1404 




1200 


1200 




210 




90 


8,2,2,3 
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TABLE 3 
Derived Parameters for L and T Dwarfs (0.9—14.5 /^m) 



Object 



Spectral Type^ 
Optical Infrared T^s (K) 



Atmospheric Parameters 
By Spectral 1^'ittini 



(cm s ^) J^ /mc 



By Evolutionary Sequences Golimowski et al. T(,g (K) 
logg (cm s~ ) 



0.1 Gyr 


3 Gyr 


10 Gyr 


1950 


2250 


2275 


1950 


2250 


2275 


1650 


1900 


1975 


1475 


1750 


1800 


1475 


1750 


1800 


1175 


1425 


1475 


1150*^ 


1375 


1425 


1150 


1425 


1500 


1150 


1425 


1500 


1150 


1425 


1500 



2MASS J1439+1929 LI LI 



2MASS J1506+1321 
2MASS J0036+1821 
2MASS J2224-0158 
2MASS J1507-1627 
2MASS J0825+2115 
DENIS 0255-4700 
SDSS J1254-0122 



2MASS J0559-1404 T5 T4.5 



L3 




L3.5 


L4±l 


L4.5 


L3.5 


L5 


L5.5 


L7.5 


L6 


L8 


L9 


T2 


T2 



2100 
2100 
1800 
1700 
1700 
1700 
1400 
1400 
1200 
1200 
1200 



5.0 
5.0 
4.5 
5.5 
4.5 
4.5 
4.5 
5.5 
5.0 
5.5 
5.5 



2 


0.582 


3 


0.418 


1 


1.000 


3 


1.000 


1 


0.999 


2 


1.000 


1 


1.000 


2 


0.988 


3 


0.831 


2 


0.105 


4 


0.999 



5.188 
5.188 

4.924 
5.378 
5.451 
5.482 
5.547 
4.799 
4.908 
4.713 



Note. — Only models with /mc > 0.1 arc listed. 
" Spectral types of the L dwarfs are from lKirkpatrick et "a!] J199gl l. lReid et ail l!2000'). Gizis ct al. ( 2000). Kirkp atrick et al.l ll2000D . IBurgasser et al.N2006all . 
and J.D. Kirkpatrick (in preparation). Spectral types of the T dwarfs are from Burgasscr ct al. (2003a, 2006a). Errors on spectral types are ±0.5 subclass 
unless otherwise noted. DENIS 0255—4700 was not included in the lGolimowski et al.| (120041) sample. We have computed its Lboi using the new parallax 
measurement of ICosta et al.l (1200611 ■ The bolometric flux was determined by integrating over its 0.6 to 14.5 /.im spectrum as described in ICushing et al.l 



(Unni). 



Model Comparisons 43 



TABLE 4 
The Properties of SDSS 1254-0122 and 2MASS 0559-1404 





SDSS 1254-0122 


2MASS 0559-1404 


Parameter 


single 


equal 


mass binary 


single 


equal 


mass binary 


log (ibol/i©) 

Tcff (K) 
logt; (cm s ^) 
Age (Gyr) 


-4.60 
1200 
4.71 
~0.3 




-4.90 

<1160 

<5.38 

<10 


-4.63 
1200 

4.80 
~0.4 




-4.93 
<1150 
<5.38 
< 10 











TABLE 5 








t 






Derived Parameters of 


THE L AND T DWARI 


7S FROM Restricted Spectral Intervals 








Object 


Spectral Type'' 






Teff(K) / log3(cms 


^) / /sod / /mC 










Optical 


Infrared 


Y 


J 


K 


L 


IRS/SL 


Near-lntrared 




2MASS 1439+1929 


LI 


LI 


1700/5.5/3/1.000 


2400/4.5/1/1.000 


2100/5.0/3/1.000 


2400/5.5/2/0.922 


2200/5.5/3/0.922 


1800/5.0/3/1.000 




2MASS 1506+1321 


L3 




1800/5.5/4/1.000 


1500/5.5/1/1.000 


1900/4.5/3/0.761 
1900/4.5/2/0.240 


2400/5.5/4/0.997 


2400/5.5/1/1.000 


1600/4.5/2/1.000 




2MASS 0036+1821 


L3.5 


L4±l 


1700/5.5/4/0.988 


1700/5.5/2/1.000 


2000/5.0/4/1.000 


2300/5. 5/nc/l. 000 


2400/5.5/1/0.863 


1700/5.5/3/1.000 




2MASS 2224-0158 


L4.5 


L3.5 


1600/5.5/3/0.996 


1500/5.0/1/1.000 


2400/5. 5/nc/0.748 
2300/5.0/nc/0.133 


1900/4.5/4/0.219 
2200/5. 5/nc/0.190 
2300/5. 5/nc/0.122 

1800/5.0/4/0.114 


2400/5.5/1/1.000 


1400/4.5/1/1.000 




2MASS 1507-1627 


L5 


L5.5 


1600/5.5/3/1.000 


1500/5.5/2/0.999 


1700/5.0/3/1.000 


2300/5. 5/nc/l. 000 


1900/5.0/3/0.627 
2000/5.0/3/0.370 


1700/5.5/3/1.000 




2MASS 0825+2115 


L7.5 


L6 


1400/5.5/2/0.549 
1500/5.0/2/0.359 


1300/4.5/1/0.984 


1500/4.5/2/0.835 
1700/5.0/3/0.166 


1800/5.0/4/0.998 


2400/5. 5/nc/O. 201 
2200/5.5/4/0.166 

2300/5. 5/nc/0.157 
1900/5.5/4/0.142 


1400/4.5/1/1.000 




DENIS 0255-4700 


L8 


L9 


1300/5.5/2/0.992 


1300/5.5/2/0.999 


1700/5.5/4/0.820 
1800/5.5/4/0.177 


1900/5. 5/nc/0.447 
1700/5.5/4/0.241 
1700/5.0/4/0.186 


1200/5.5/1/0.465 
1500/4.5/2/0.218 
1200/5.0/1/0.167 


1500/4.5/2/0.572 
1400/4.5/2/0.428 




SDSS 1254-0122 


T2 


T2 


1000/4.5/2/1.000 


1200/4.5/3/1.000 


1500/5.5/4/1.000 


1500/5.5/3/0.918 


1400/5. 5/nc/O. 910 


1400/4.5/4/1.000 




2MASS 0559-1404 


T4.5 


T4.5 


1200/5.5/nc/l.OOO 


1100/4.5/4/1.000 


1100/4.5/2/1.000 


1100/4.5/nc/0.469 
1000/4.5/3/0.366 


1200/5.5/4/0.982 


1100/5.0/4/1.000 


o 


Note. — Only mode 


i\s with /f 
L dwarfs 


AC > 0.1 are listed. 










ind J.D. Kirkpatrick 




^ Spectral types of the 


are fromlKirkpatrick et al.' fl999"), R 


cid ct al. f2000l)JGizii 


5 et al.l (l200CH.IKirkpatrick et all (I2000D, IBureasser et al.l 12006311, e 


M 


(in preparation). Spectral types 


of the T dwarfs 


are from.Bureasscr ct al., (,2003a„ .2006a). Errors on spectral tvpef 


1 are ±0.5 subclass unless otherwise noted. 
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